Spectroscopy of the Potential Profile in a Ballistic Quantum Constriction 
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We present a theory for the nonlinear current-voltage char- 
acteristics of a ballistic quantum constriction. Nonlinear fea- 
tures first develop because of above-barrier reflection from 
the potential profile, created by impurities in the vicinity of 
the constriction. The nonlinearity appears on a small voltage 
scale and makes it possible to determine distances between 
impurities as well as the magnitude of the impurity poten- 
tials. 
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In recent experiments on microconstrictions, electro- 
statically created by a split-gate technique in the 2D 
electron gas of a GaAs-AlGaAs heterostructure, several 
groups HH have reported small-scale nonlinear features 
of the system conductance when measured as function of 
the driving voltage, V. These anomalies, which occurred 
on the first fundamental 2e 2 /h- conductance plateau, still 
lack a satisfactory explanation. In this Letter we show 
that the small-scale nonlinearity in the conductance is de- 
cisively determined by the shape of the effective potential 
barrier inside the quantum ballistic constriction (QBC). 
The latter is defined by the geometry of the microcon- 
striction and by impurities in its vicinity. The effective 
potential profile leads to an 'above-barrier' reflection ef- 
fect ||, which gives rise to the nonlinearity in the QBC 
under consideration. In short, the characteristic scale of 
the nonlinear features is determined by the characteristic 
inter-impurity distances. 

Our starting point is the well known adiabatic model 
H of a QBC shown in Fig. |]. Numerical calcula- 
tions confirm that the effect of charged impurities 
in the vicinity of the microconstriction is to smoothly 
modulate the electrostatic potential inside the constric- 
tion, Vi mp {x,y). This modulation can be thought of 
as a simultaneous modulation of the microconstriction 
geometry and a smooth variation of the effective one- 
dimensional potential for longitudinal transport along 
the QBC axis. Such a picture is valid as long as the 
condition \dVi mp /dy\ d(0) <C AE± is satisfied. Here 
AEj_ ~ n 2 h 2 /2m*d{0) 2 is a characteristic intermode 
level spacing in the adiabatic quantum constriction. 

Assuming the constriction to be long and smooth 
enough, the criterion d'(x) ~ d(0)/l <C 1 is fullfilled (I 
is the length of the contact), and we are in the adiabatic 
transport regime. Under these conditions — to zero order 
in corrections to adiabatic transport — the Schrodinger 
equation for the longitudinal part of the wave function 
in the ra-th transverse mode is 



2m 



-4>n(x) + W n (x)lp n (x) = Elp n (x). 



(1) 



Here W n (x) = Vi mp (x, 0) + E n ±(x) corresponds to an ef- 
fective gate-generated potential profile modulated by the 
impurity potential both directly (Vi mp (x, 0)) and through 
the modulation of the constriction width, which depends 
on the transverse component of the impurity potential. 

We consider an electron in the n— th propagating 
mode, and assume that the electron energy is everywhere 
much higher than the effective potential, E ^> W n (x). In 
this case above-barrier scattering will lead to a reflection 
coefficient of the form |Hl 



Rn(E) 



m*/(h 2 k(E)) / W n {x)e 2lk ^ x <k 



, (2) 



where k(E) = y/2m*E/h. 

The current through the contact is given by 

I(V) = (2e/h) J dE [n F (E - m) 

N 



(3) 



Here N is a number of propagating modes, jii — Ep — 
j3eV and fi2 = Ep + (1 — (3)eV are the chemical potentials 
of the left- and right electron reservoirs, and < (3 < 1 is 
factor that describes the asymmetry of the voltage drop 
along the constriction 

The most spectacular situation occurs at zero temper- 
ature, when a simple expression for the differential con- 
ductance of the system follows from Eq. (|^). One finds 
that 



G(V) = N(2e 2 /h)-5G N (V), 



(4) 



where 



N 



5G N (V) = (2e 2 /h) Y, [PRn(E F - 0eV) 

71=1 

+ {l-p)R n (E F + (l-(3)eV)}. 



(5) 



Clearly, the reflection coefficient (g) is determined by 
two factors; first by the contact geometry, and second by 
the impurity potential. The geometric contribution to 
W n (x), i.e. E n j_(x), is a sum of a smooth function and 
a small-scale oscillatory term. The former is due to the 
gate potential, while the latter is induced by the trans- 
verse part of the impurity potential. It is well known 
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that the smooth part of E n ± (x) contributes an exponen- 
tially small term to R n 0). Nevertheless, as the number 
of propagating modes increases, the total amplitude due 
to this contribution grows as iV 5 (as E n ± « n 2 ). There- 
fore, if only the number of open modes becomes large 
enough, the relative amplitude of the oscillatory compo- 
nent of the reflection coefficient is weakened. The best 
conditions for doing 'potential profile spectroscopy' are 
then at hand when the number of propagating modes is 
as small as possible. 

From now on we take N — 1. We can thus neglect the 
contribution to R n from the adiabatic part of the trans- 
verse energy that enter W n (x) . Evidently, the other con- 
tributions to this function, coming from the irregularities 
of the constriction geometry and from the potential pro- 
file inside the constriction, can not be separated. This 
corresponds to the physical situation, since both these 
contributions have the same origin: the random poten- 
tial of the adjacent impurities. We denote this part of the 
potential W n {x) — which is independent of mode index 
n — by Vim P (x). 

Note that by considering only a single propagating 
mode we automatically exclude opening/closing of con- 
ducting modes Jl],|9|,[l0| as well as coherent mode mixing 
p| as possible mechanisms for the nonlinear features 
discussed. 

The effective potential Vi mp (x) is written in the form 



Kmp(^) ^ ^ U a (\x X a \ 



(6) 



where U a denotes the contribution to the self-consistent 
effective potential from the a-th impurity, situated at 
the point x a . The reflection coefficient then acquires the 
simple form 



R(E) = (m* /Ti 2 k{E))' 



Y^e 2lk{E)x -u a {E) 



(7) 



u a {E) 



dxe 2lk{E)x U a {\x\) 



where u a (E) is the Fourier transform of the effective im- 
purity potential. 

It is evident from (^) that the reflection coeffitient 
is sensitive to the interimpurity distances. It follows 
from the uncertainty condition for the direct and in- 
verse Fourier transform that the accuracy Ax with which 
these distances can be determined from (0) is limited by 
AxAk > 1/4-7T. In our case this reduces to 



Ax >~ A FJ B F /(47r 2 |eV r |). 



(8) 



The small numerical factor \/Att 2 allows for relatively 
small driving voltages to define the interimpurity dis- 
tances with accuracy of order Xp. 



It is easy to estimate from Eq.(0) the order of magni- 
tude of the relative conductance variation due to above- 
barrier scattering at zero temperature: 



5G_ 
~G 



J 1 imp Uimp 
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Because of the smoothness of the effective impurity 
potential the quantity UmpUimp is the maximal value of 
its Fourier transform, limp being its effective length, and 
Uimp its magnitude. The characteristic scales for these 
quantitites are given by different authors in a rather wide 
interval. Following ||, we use Ui mp /Ep ~ 0.3; li mp /Xp 
can be estimated to be in the range 0.1 to 0.01. That 
is, the magnitude of the nonlinear contribution to the 
conductance at zero temperature is 1 — 10 _2 %. 

At finite temperatures the conductance oscillations will 
be washed out. Thus in order to observe the effect de- 
scribed it is necessary that UbT <C eAV, where AV is 
the period of oscillation with bias voltage. A realistic es- 
timate is that the effect can be observed if T < lOOmK. 

Numerical calculations of the conductance vs. bias 
voltage are presented in Figs. ^- ||. We use Coulomb 
impurity potentials, 



U a (\x - x a \) = Z a e 2 / «V (x - x a ) 2 + h? a . 



(10) 



The second term under the square root determines the 
lateral and vertical displacement of the a-th impurity 
from the axis of the constriction. The dielectric constant 
for GaAs is K ~ 13. An asymmetric voltage drop {(3 ^ 
0.5) leads to an asymmetric response SG(V). Evidently 
the larger the inter-impurity distance, the more distinct is 
the picture of oscillations, and the larger is the accuracy 
with which we can determine the potential profile (see 
Eq.®). 

The average value of the nonlinear term in the conduc- 
tance is described by the averaged value of the function 
R(E) ■ E (see Eq. (||)). Making use of the Parseval in- 
equality, we obtain the result 



1 r°° 

(R(E)E) = — / dkR(E)E 



Ahr J- 



dx 



a 

(ii) 



In this way we find the average magnitude of the impurity 
potential along the axis of the constriction. 

Since the Fourier transform of R(k(E))E directly gives 
the set of interimpurity distances, Ax aa > , we can apply 
a filtration procedure on the corresponding frequencies, 
that is find 



1 f°° 

{RE). = — / dkR(k)E(k)exp(-2ikAx aa > 
2vr J^^ 



(12) 
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Hence we can analyze the lineshape of every peak in the 
trace of 6G(V), initiated by each distinguishable impu- 
rity. In the general case, such an analysis allows us to 
determine the quantities \h a — h a > \ / \h a + h a i | , thus pro- 
viding information about the distances between the im- 
purities and the constriction plane. 



Using the model potentials (10) we find for the line 
shape an ex- 

pression (RE) Ax / = (m*/h 2 )(e i Z a Z a ,/K 2 \h a + h a >\) ■ 
K{\h a -h a >\/\ha+h a >\). The function K(£) is the 
complete elliptic integral of the first kind. 

The nonlinear contribution 6G(V) to the conductance 
allows us to determine the spectral characteristics of the 
potential profile created by impurities in the constric- 
tion. For real impurity potentials U a (k) are smooth 
functions of the variable k and therefore of the driv- 
ing voltage V. Because of this fact we can apply the 
Fourier transformation technique to Eq. (^) and from the 
positions of the spectral peaks determine the distances 
Ax aa > = \x a — x a '\ between the impurities. Evidently, 
these distances can be measured only if they are larger 
than Xp. If on the contrary Ax aa > < Xf we can regard 
the corresponding impurity cluster as a single scatterer. 
Anyhow, it follows from (||) that they can be resolved 
only at large enough voltages V, i.e., in the regime where 
the nonlinear effects are developed. A Fourier analy- 
sis of G(V) as well permits us to find the quantities 
H E Q \ U <*\ 2 and J2 afa , U a U*, cos(Arg([/ Q [/*,)), which 
characterize the magnitudes of the impurity potentials. 

In Figs||-|] different types of conductance- voltage 
curves caused by above-barrier reflection in the constric- 
tion are shown. We plot there the nonlinear contribu- 
tion 5G to the differential conductance(in arbitrary units) 
vs. the reduced driving voltage v = V/Ep. The scale 
of v G [—0.5; 0.5] is consistent with measured values of 
Ep ~ lOmeV and distance between the first and second 
subbands in a QBC, AE± ~ 5meV [1]. The potential 
drop is assumed to be completely asymmetric (f3 = 1). 

In Figj^ the next simplest case is shown, with a single 
impurity. Evidently the nonlinear correction is a mono- 
tonic function of voltage. Nonmonotonitic behavior ap- 
pears when two impurities are present (Fig j3j) . Fig.^ 
demonstrates the case of three impurities. It is easy to 
see that the structure of the curves becomes richer as the 
number of impurities and/or average distance between 
them grows. 

In conclusion, we have demonstrated that above- 
barrier scattering of electrons by the potential profile in 
quantum ballistic constrictions leads to small-scale non- 
linearities in the current- voltage characteristics. The ef- 
fect is best observed when only one conducting channel 
is open, and at temperatures less than ~ 100 mK. It 
can be used to determine the spatial arrangement of the 
impurities in the constriction. 
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FIG. 1. Quantum ballistic constriction, (a) Clean adia- 
batic channel; d(x) is a smooth function of longitudinal co- 
ordinate, (b) Effective potential profile inside the channel in 
presence of impurities. 



FIG. 2. Nonlinear contribution to the conductance of the 
constriction with a single impurity. 
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FIG. 3. Nonlinear contribution to the conductance of the 
constriction with two impurities. The oscillations arise due 
to Fabry - Perot type interference of electron wave between 
them (the distance between the impurities is 11\ F ). 

FIG. 4. Nonlinear contribution to the conductance calcu- 
lated for 3 impurities in the constriction. The average dis- 
tance between the impurities equals to 23 Af . 
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